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Abstract 

This work uses Lorentz-signature in-in perturbation theory to analyze the late-time 
behavior of correlators in time-dependent interacting massive scalar field theory in de 
Sitter space. We study a scenario recently considered by Krotov and Polyakov in which 
the coupling g turns on smoothly at finite time, starting from (7 = in the far past 
where the state is taken to be the (free) Bunch-Davies vacuum. Our main result is that 
the resulting correlators (which we compute at the one-loop level) approach those of the 
interacting Hartle-Hawking state at late times. We argue that similar results should 
hold for other physically-motivated choices of initial conditions. This behavior is to be 
expected from recent quantum "no hair" theorems for interacting massive scalar field 
theory in de Sitter space which established similar results to all orders in perturbation 
theory for a dense set of states in the Hilbert space. Our current work i) indicates 
that physically motivated initial conditions lie in this dense set, ii) provides a Lorentz- 
signature counter-part to the Euclidean techniques used to prove such theorems, and 
iii) provides an explicit example of the relevant renormalization techniques. 
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1 Introduction 

Quantum field theories (QFTs) in de Sitter space are of interest for many reasons. Some 
of these include the increasingly precise measurements of the cosmic microwave background 
(CMB) [ ] which have prompted many to study predictions of the CMB spectrum beyond 
the Born approximation (see, e.g. [2, 3, 4, 5, 6, 7, 8, 9, 10]). Others include a growing 
interest in understanding local measurements in eternal infiation (e.g., [11, 12, 13]), as well 
as a renewed interest in approaches to de Sitter quantum gravity [14, 15, 16, 17] inspired 
by dS/CFT [18, 19]. In addition, the fact that de Sitter (dS) is a maximally symmetric 
(and thus relatively simple) example of a spacetime where horizons limit the observations of 
freely-falling observers makes QFTs on dS of interest in their own right. 

One of the chief concerns with QFTs in de Sitter has been their infrared stability (see e.g. 
[20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40]). Two recent 
papers [38, 39] have established that at least massive scalar field quantum field theories are 
infra-red stable (at all orders of perturbation theory) in a particular sense. In order to state 
their results precisely, we first introduce the (interacting) Hartle-Hawking state |HH) [11] 
defined by analytically continuing all correlation functions from Euclidean signature. Next, 
consider a normalized state constructed by the application of smeared field operators on 
|HH): 

1^)=/ /(2/i,...,Z/„)<^(l/i)-->(2/„)|HH), (1.1) 
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with f{yi, . . . ,yn) a smooth smearing function of compact support. By the Reeh-Schheder 
theorem of curved spacetimes [12] the set of states of the form (1.1) is dense in the Hilbert 
space. The works [38, 39] show that, at all orders of perturbation theory, the correlation func- 
tions of |\E') reduce to those of the Hartle-Hawking state when evaluated in the asymptotic 
future/past of de Sitter space: 

{(j){Xi) ■ ■ ■ (j){Xn))^ -> {(j){Xi) ■ ■ ■ . (1.2) 

In particular, de Sitter invariance of |HH) means that all one-point functions approach 
constants (whether the associated operators are elementary or composite). Ref [■>>;] calls the 
result (1.2) a 'quantum cosmic no hair theorem,' while in the language of [■'> 1] one says that 
|HH) is an attractor state for local correlators. 

Although stated as a result concerning QFTs in exact de Sitter space, the no-hair theorem 
just described may also be usefully applied to more interesting scenarios. Since it relies only 
on the asymptotic behavior, it should be valid in the asymptotic region of any asymptotically- 
de Sitter spacetime, or within the causal patch of an observer who finds herself in a locally 
de Sitter spacetime. One expects physically relevant states to take the form (1.1) within the 
de Sitter region so that the theorem applies. 

Our main purpose here is to provide evidence that this is indeed the case by studying 
(at the one-loop level) a particular scenario recently discussed by Krotov and Polyakov [10]. 
In this scenario, the spacetime is again exact de Sitter but the theory is time dependent. 
The particular model involves a cubic interaction g{x)(j)^{x) with time-dependent coupling 
g{x), taken to vanish in the far past and to approach some constant gj in the far future. 
The state of the system is taken to be the free Bunch-Davies vacuum in the region where 
g{x) = and we take the coupling to turn on at some fixed time. We explicitly compute 
the 0{g'^) (one- loop) corrections to the 2-point function in this model and verify that they 
approach those of the de Sitter Hartle-Hawking state in the far future. We work entirely in 
Lorentz signature, in part to counter concerns [26, 32, 40] about the Euclidean techniques 
used in [38, 39]. In addition, we note that the current work provides an an explicit example 
of the renormalization techniques used in [ >' )] which combine Pauli-Villars regularization 
with Mellin-Barnes representations. Our techniques also apply to other examples where the 
coupling does not depend on time but in which the spacetime is de Sitter only after some 
finite time. 

Although [ ''■] concluded that de Sitter QFT is "unstable," it is useful to point out that 
our technical results are completely consistent with those of [l!;]. As stated below their 
equation (17), for fixed g{x) their approximations are not valid for correlators computed at 
late times. Instead, [ 10] focused on correlators defined at some fixed time in the limit where 
g{x) turns on at very early times. The divergence they find is in fact to be expected from 
the results of [38, 39], which suggest that correlators in well-behaved states approach those 
of |HH) in the far past. Using the free vacuum when g{x) = and taking g{x) to turn on 
very early ensures that correlators at such early times differ significantly from those in |HH). 
As a result, one already expects from [38, 39] that the state is not well-behaved under the 
limit taken in [40]. 
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We emphasize that the present work considers only perturbative effects in massive the- 
ories. Non-perturbative effects can yield qualitatively different behavior (though see [13]), 
and including massless fields (whether scalar or tensor^) would raise new issues. As a result, 
our analysis does not directly address the much- discussed possibility of novel infrared effects 
in de Sitter quantum gravity (e.g., [46, 47, 48, 49, 50, 51, 52, 53, \), 10]). Nevertheless, a 
detailed study of theories on a fixed de Sitter background helps to disambiguate those effects 
which are truly quantum gravitational from those that are generic for quantum fields in de 
Sitter. 

We begin in §2 with a brief review of linear quantum field theory on de Sitter. We then 
analyze time-dependant couplings in global de Sitter in §3. Section 3.1 studies a simple 
theory with a time-dependant 0^(x) interaction (i.e., a time-dependant mass perturbation). 
We study the more complicated model of a time-dependant (/)^(x) interaction in §3.2, relying 
at times upon the results derived in the 0^(a;) model. Some further calculational details are 
presented in appendices. We provide a concluding discussion in §4. 



2 Preliminaries 

We begin by reviewing some basic features of free quantum fields in de Sitter space. Recall 
that global de Sitter may be described by the metric 



^ dr]^ + {l + ri^)dnl_. 



(2.1) 



where i is the de Sitter radius, r/ is a time coordinate with range — oo < r] < +oo, and dfljj_i 
is the metric of the unit sphere S^'^. The coordinate i] is related to the more familiar global 
de Sitter coordinate t (for which gu = —1) via r] = sinh(t/£). In these coordinates the 
volume element is ^^/—g{x)d^x = {1 + ri'^)^^~'^^^'^dri d^ln-iix) with x a unit vector in ]R^~^ 
parameterizing S^^^. The Euclidean section of de Sitter is the Euclidean sphere with 
radius i. 

Free massive scalar field theories on de Sitter have been well-understood for decades (e.g., 
[54, 55]). Such theories may be described by the classical Lagrangian density 

1 

C = - V^0V^0 + — 0^ (2.2) 

from which it follows that the classical equation of motion is the Klein-Gordon equation; 
correspondingly, in the quantum theory the Schwinger- Dyson equations are 

(□, - M^) (0(a;i) ■ ■ ■ 0(x,) ■ ■ ■ 0(x^))^ = 0. (2.3) 

The free theory admits a unique Hadamard de Sitter-invariant state known as the Bunch- 
Davies, Euclidean, or (free) Hartle-Hawking state, which we denote by |0). The latter two 



^The 2-point functions of Maxwell fields in dS are known to behave like those of massive scalars [44, 4-1] 
and so provide no new subtleties. 
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names come from the fact that the correlation functions of this state may be defined by the 
analytic continuation from the Euclidean section. We denote the time-ordered, anti-time- 
ordered, and Wightman 2-point functions of this state by 



Wa(Xi,X2] 



(T0^(Xi)0<^(x2))o , 

(T0^(Xi)0<^(x2))q 

{(j)a{Xl)(j)a{X2))o. 



(2.4) 



In these expressions we have introduced a label a to keep track of the bare mass M. We 
define a by 



D - 1 



+ 



4 



M 



2n2 



1/2 



(2.5) 



from which it follows that M'^P' = —a{a + D — 1). 

At times it will be convenient to expand the scalar Green's functions (2.4) in Klein- 
Gordon modes 0^|;(x). These modes are orthonormal with respect to the Klein-Gordon 
inner product: 



fjD-l, 



tL' v^o-m/ kg •" 



(l + ^2)(D-i)/2 f dnn-,{x)n^ U £(r/,f)V,e,^(r/,: 
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LM- 

(2.6) 



Here is the future-directed normal vector [n^ = (1 + ff y/'^d'^ / tj to an 77 = const, surface 
and aV^B :^ AV - BV ,A. The KMn-Gorion .ode. .nay be written expUcitly a. 



with Yj^{x) spherical harmonics on S^~^ and Urjiiv) given by 

(L+(D-2)/2)/2 



2^-{D-2)/A 



1 — irj 
1 + ii] 



where F^l^t]) is a Gauss hypergeometric function 

D 



Fadv) = 2F1 



D ^ D l-iri 

— ; ^ H — ; 

2 ' 2 ' 2 



and the normalization coefficient is 

7V„ 



T[L-a,L + a + D-l] 



1/2 



" r(L + f) 

Using these modes we may expand, e.g., the Wightman function 

W^{X,,X2) = £2-^5^0^£(Xi)0;^(x2) 



(2.7) 



(2.^ 



(2.9) 



(2.10) 



D-2\ 00 



ff2-D 



4t^d/2 



Y^{2L + D- 2Mvi)<L{V2)Ci''-'^^\x, ■ X2). (2.11) 



L=0 
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Here C 



{D-2)/2, 



z) is a Gegenbauer polynomial and once again xi, X2 are unit vectors in 



parameterizing the ^. To obtain the last equality we sum over angular momenta via the 
useful identity 



(2L + D-2)Cf-')/'(fi-X2), 



(2.12) 



It will be useful to note some qualitative features of the of the Klein-Gordon mode 
function Ucriiv)- First, as one might expect from the fact that the volume of the S^~^ is 
smallest at 77 = 0, the mode functions are bounded by their values at rj = 0: 



KL(r])r<KL(0)r = 2-(2^+^-i)7r 



T[L-a,L 



a 



D-l] 



(r [i 



+L-0- I+L+0-+D-I 



])• 



2L 



(l + 0(L-i)) 



2 ' 2 
when L 1, a fixed 



(2.13) 



(see eq. (56) of [jG]). From this we see that \u„l{vi)\ may be bounded by a function of L 
that decreases as L — )• cxd. 

Second, the expansion of the universe and the ensuing growth of the physical wavelength 
at fixed L suggest that at large \rj\ all modes behave like the L = mode. Indeed, as derived 
in detail in Appendix A, the following asymptotic expansion for Uaiiji) is valid for large |?7| 
when cr, L, and D are held fixed: 



2(T+{D-2)/2 



L+^,2a + D-l 
L + (T + D-l,a+^ 



exp 



TT 



^-(L + cT + D-2)\ {r^r 



1 + 



[a ^ -{a + D -1)), for |r/| > 1, jr^] > (L - (t). 



Now, while (2.14) gives the correct asymptotic behavior for a given mode at large \ri\ 
(with all other parameters fixed), it does not correctly reproduce the behavior of the mode 
function at some arbitrarily large-but-finite |?7| as L — ?■ 00. To understand the behavior of 
the mode function in this regime, we instead use the WKB approximation valid when 

> 1. (2.15) 
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The WKB approximation is derived in appendix B and is given by 

1 



g±iT(T,) 



with T{ri) satisfying 



d 

dr] 



T(r^) 



fiv) 



.1/2 



(2.16) 



(2.17) 



The key feature of this expression is that T{ri) is large in the regime of validity, so (2.16) is 
a highly oscillatory function of rj. 
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•X2 



n 



Figure 1: The Penrose diagram of de Sitter. We consider the state defined by the Bunch- 
Davies vacuum on a Cauchy surface S. The interaction turns on across the region TZ via a 
smooth couphng function g{x) such that g{x) = gf in the causal future of 71 and g{x) = 
in the causal past of TZ. We compute the time-ordered 2-point function (T0(xi)0(x2))^ of 
two points in the distant future. 

Let us now turn the discussion to interacting theories. The Hartle-Hawking state |HH) is 
constructed by analytically continuing all correlation functions from the Euclidean section. 
We denote the Hartle-Hawking state constructed perturb at ively in an interacting theory by 
|HH), and reserve |0) to denote the Hartle-Hawking state of the free theory. The state |HH) 
has been studied in detail for massive scalar field theories in [30, 34, 38, 39, 57] by performing 
the relevant analytic continuations. However, in the current work we perform all calculations 
explicitly in Lorentz-signature using standard Schwinger-Keldysh (a.k.a "in-in", "real-time", 
"closed time path") perturbation theory (for original works see [58, 59]; for more tractable 
introductions see [60, 61, 62] and the appendix of [j]). 

3 Time-dependent couplings in de Sitter 

Consider a massive scalar field on global de Sitter with a time-dependent self-interaction. In 
particular, let the self-interaction vanish in the asymptotic past but turn on smoothly across 
a spacetime region 71. We require the coupling function g{x) to satisfy 



where as usual denotes the causal past and future of a set and A/B denotes the set of 
points in A that do not lie in B. We sketch the scenario in Fig. 1. 

We wish to compute correlation functions with respect to the state |\E') which coincides 
with the free Bunch-Davies vacuum |0) in the past region J^{7V)/7Z. For infinitesimal 




for X G J-(7^)/7^ 
for X e J+(7^)/7^ 



(3.1) 
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Xi • X2 Xi • • X2 

(a) (b) 

Figure 2: Corrections to the time-ordered 2-point function (T0(xi)0(a:2))^ due to the in- 
teraction — ^(yf(x)0^(x). The 0{g{x)) correction is depicted in Fig. (a) while the 0{g'^{x)) 
correction is depicted in Fig. (b). It is convenient for computation to label each leg of the 
diagram by a distinct mass parameter cTj. 

coupling g{x) the time-ordered 2-point function with respect to |\E') can be expanded as 

00 

(T</)(xi)0(x2))^ = {T<p{x,)<f){x2)), + J2 {T<j>{x,)<j){x2)}P , (3.2) 

n=l 

where (T0(a;i)0(x2))^'' is of 0{g'^). We choose a surface S in the past of TZ as our initial 
Cauchy surface - see Fig. 1. Given the above choice of state, the appropriate Green's 
functions to use in the Schwinger-Keldysh formalism are those of the Bunch-Davies vacuum 
(2.4). 

The main result of this section is to show that, when it is evaluated at xi, X2 in the far 
future of 7^, the 2-point function (3.2) reduces to that of the Hartle-Hawking state of the 
analogous theory with constant coupling g^. We consider both quadratic and cubic couplings 
{g{x)(j)'^{x) and g{x)(j)^{x)) below, showing in each case that the perturbative corrections 
(T0(a;i)0(x2))|^'* approach those of the Hartle-Hawking state at late times up to order n = 2. 
Although it may be of less physical interest, the analysis of the quadratic coupling in section 
2 will help to greatly simplify our 1-loop treatment of the cubic coupling in section 3. 

3.1 Example: interaction 

In our first example we consider the simple quadratic interaction term 

Ant[0] = -^0'(X). (3.3) 

Treating this term perturbatively yields only tree-level diagrams; no regularization or renor- 
malization is needed. We compute both the 0{g) and the 0{g'^) corrections to the two-point 
function of this model below. As we will see in section 3.2, the 0{g^) correction from (3.3) is 
closely related to the 1-loop correction from a 0^ interaction. As a result, the results below 
will greatly simplify the manipulations in section 3.2. 
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3.1.1 The 0{g) correction 

The 0{g) correction is depicted in the Feynman diagram shown in Fig. 2 (a) and is given by 
the expression 

{T(f)^^{xi)(p^^{x2))^^^ = i / g{y) [G„,{y,xi)G^^{y,X2) - Wa^{y,xx)W^^{y,X2)\ ■ (3.4) 

We denote by ^y - ■ ■ an integral over the future of S. For the moment it is convenient to let 
each Green's function have a distinct mass; we will take the limit of equal masses later. 
Consider the first term in (3.4): 



Tx{xx,X2) :=% / g{y)G„^{y,xx)G„^{y,X'^. 

J V 



(3.5) 



Making use of the Green's functions' equations of motion 

(□, - M2)G,(x, y) = (□, - M')G^{x, y) = i5{x, y), 
(□,. - M')W^{x, y) = (□, - M^)W„{x, y) = 0, 

we may usefully re-write (3.5) as 

i 



(3.6) 



Ti{xi,X2) 



9iy) 



[ayG^^{y,xi))Ga.2iy,X2) 



Mf-MlJy 

-GaAy,xi){ayGa2iy,X2)) - i5{xi,y)G^^{y,X2) + i5{x2,y)G^^{y,xi) 
[g{xi)Ga2 {xi, X2) - g{x2)Ga^ (xi, X2)] 



Ml - Mi 

i 



Mf - M| 



9{y) 



[^yG^^{y,xi))G„.,{y,X2) - G^^{y,xi){UyG^^{y,X2)) 



(3.7) 



Making the simple re-arrangement 

g{.y){.^G„-,{y,xi))G„^{y,X2) = W[g{y){V f^G^^iy, xi))G^^{y, X2) 

-iV^g{y)){V,G^M^i))G.M^2) 
-g{y){V^G^,{y,x,)){V^G^,{y,X2)), 



(3. 



we obtain 



Ti(xi,X2) = ^2 _ J^2 [9iXl)Ga.,{Xi,X2) - fi'(x2)G^i(Xi,X2)] 



M!-Ml,y 



X2 



Ml - Ml Jy 



(V^^(l/)) [G,,(i/,xi)V^G,,(y, 



X2) 



(3.9) 
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One can perform the same manipulations for the second term in (3.4). The only difference 
is that the Wightman function satisfies the homogeneous equation of motion, so there are 
no analogs of the terms on the top line of (3.9). All together we obtain the expression 



{T(f)^^{xi)(f)^^{x2)y^ 

i 



(1) 



[g{xi)G„^{xi,X2) - g{x2)G„^{xi,x2)] 



i 



Ml - M| 

g{y) (^Ga^{y, xi)v ^G^^^y, X2) -Wa^{y,xi 



X2, 



Mt - Mi jy 



iV^g{y)) \G^,{y,x,)V^G^,{y,X2)-W,,{y,x,)'^^W^,{y,X2)] . (3.10) 



In the second line of (3.10) there is an integral of a total derivative. By Stokes' theorem 
this integral can be expressed as an integral over the boundary of the region to the future 
of S. This boundary is simply the union of S and future infinity /+. Now, the combination 
of Green's functions in the integrand is such that the integrand has support only on the 
union of the past light cones of xi and X2, so the integral over /"•" vanishes. Furthermore, the 
coupling function g{y) vanishes on S, so the integral over S vanishes as well. We conclude 
that the integral in the second line of (3.10) is identically zero: 



(T0^,(x 1)0^2 (^2))i'^ 



r2 nj2 [9iXl)G^2{Xi,X2) - gix2)G^,{Xi, X2)] 



Ml - M| 



Ml - Mi 

(V^5'(?/)) G^^{y,Xi)V f,G^.,{y,X2) - W^^{y,Xi 



.(3.11) 



We are interested in xi,X2 in the future region J^{7V)/TZ. Since the gradient W^g has 
support only in TZ, we need not allow y in (3.11) to coincide with Xi, X2 or to lie in the future 
of either point. We may therefore replace Gc^iGa^ by appropriate Wightman functions in 
the integral and write the second line of (3.11) in the form 



^2 CTl 0-2(^11 ^2) 



Ml - M| Jy 



iV'^giy)) \w^Ax^,y)'^ ^W„,{x2,y) - W^.iy^xi)"^ ^W,,{y, 



X2) 



Ml - Mi 



lm\ [ {V^g{y))W^^{y,x,)'^,W^,{y, 



X2, 



(3.12) 



To keep the computation simple we choose g{y) to be a smooth function of the time coordi- 
nate ?7 alone, i.e., 

V^g{y) = -r\l + r^')g'{r^)5^^. (3.13) 
By expanding the Wightman functions in Klein-Gordon modes as in (2.11) we compute 



T, 



21 



D-2 



2cti(72 12;i, X2) 



Ml -Mi 



ImU rfr/(l + r72)(D-i)/2^/^^) f dno-iW^M ^i)'^W^2iy, 



f/2-D 



27r^/2 



L=0 



Im 5^ (2L + D - 2)x.,.2(^) K,Liti)u:^Lit2)G^ 



(D-2)/2, 



,X2, 



Xi ■ X2 
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with 



Xo-lO-2 



{Ml 



(i?7 (1 + ?7 



2n(D-1)/2 



(3.15) 



where 9„ denotes a derivative along the unit (future-pointing) timehke normal to the surface 
Tj = const. The integral in (3.15) is guaranteed to be finite for any pair u^^Lirj), u^^Lirj) as 
the harmonics are bounded and g'{ri) is smooth with compact support. It is also finite in 
the limit M| — )• M^, as can be seen by using I'Hopital's rule. 

It remains, however, to determine the convergence of the sum over L in (3.14). This is 
governed by the behavior of Xaia2i.L) large L ^ 1. Since g'{ri) has compact support it 
must vanish for \rj\ greater than some |?7|max, and for L(L — D + 2) ^ 1 + r/^^^ we may use 
the WKB approximation (2.16). Since the phase T(?7) in (2.16) is highly oscillatory at large 
L, we expect Xcricr2(-^) to decrease rapidly with L. To verify that this is the case, suppose 
that in fact L{L — D + 2) ^ M^£^(l + rj'^ax) so that we may also expand the phases Ta{v) 
(where we have added the label a to indicated the dependence on mass) as 



= ^ (L{L -D + 2))(^-")/2 



(3.16) 



n>0 



Noting that Tq.o- is independent of a, we now introduce a new time coordinate fj defined by i) 
is a smooth strictly increasing function of r/, ii) fj = To,cr in the region where g'{r]) ^ 0, and 
iii) fj = 1] at large \r]\. Then for large L (3.15) becomes {L~^ times) the Fourier transform of 
a smooth L-independent function of fj and, as a result, decays faster than any power law in 
L. (Here we use the fact that the sub-leading terms in the WKB expansion correct (2.16) by 
multiplying (2.16) by functions which become essentially constant at large L for \r]\ < rjmax-) 
For later use we note that for L ^ \rjmax\Ml we have shown that |xo-i,o-2(-^)| < 1/L"' and, as 
a result, that for all L and n we have a bound 



\xUL)\<Cn 



my 



(3.17) 



for appropriate C„ determined by only g{x) (and in particular, for which the C„ do not 
depend on M, L). 

It follows that the sum over L in (3.14) is absolutely convergent and yields a finite result, 
even when xi = X2- To see this, we bound the harmonics |M^2,(a;i)M*^(r/2)| < |w<jl(0)P as in 
(2.13), and we bound the Gegenbauer polynomial by it's value at coincidence: 



C 



(D-2)/2 



(1) 



L + D-2 
L + l,D-2 



(3.18) 



(see eq. (44) of [ -i;]). For L ^ 1 this behaves like L^^^. From these bounds it follows that 
the summand in (3.14) may be bounded at large L by L^~^XcTia2{L). Since Xo-io-2(-^) decays 
faster than any polynomial as L — oo, the sum is absolutely convergent. 

Furthermore, we can show that when |?7i|, \r]2\ S> \v\ma.x the expression T2aa{xi, X2) decays 
like {rjirj2Y . For r/i, 772 in this regime, let us choose some Lcut such that min(|?7i|, |?72|) 3> 
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4L^y^ ^ l^lmax- We then split the sum over L in (3.14) into two parts: One part T2_<(xi,X2) 
is the finite series containing the terms below Lcut- The other part T2,>(xi, X2) is the infinite 
series containing the terms with L > Lcut- For each term in T2_<(xi, 0:2) we may approximate 
W^^iVi) KtL^V^) by the asymptotic expansion (2.14). It follows that each term in the 
series T2_<(a;i, 0:2) decays like {rjirj2Y , and so T2_<(a;i, 0:2) does as well. 

On the other hand, we can bound the contribution of the infinite series T2^>(a;i, X2) by 

T2>,(a:l,X2)<^^-^-^^ (2L + D-2)|x.,.,(L)||n.,(0)rcf-^)/^(l) . 

^ L = Lr ^ ^ ^ 



(1) (3.19) 



This series converges absolutely and (due to the rapid decay of Xo-i<t2(-^)) decreases faster 
than any power law as Lcut is increased. Thus, as \rji^2\ — )■ C)0 we may increase Lcut (taking 
Lcut to be, e.g., any geometric mean of Ir/lmax and min(|r7i|, \rj2\)) and the contribution due 
to T2^>(xi,X2) becomes negligible. We conclude that the full expression T2a,a{xi,X2) decays 
like {rjirj2Y . 

For later use, let us also consider T2o-,/^(a;i, 0:2) for all Re/x < Re a. Of course, the same 
argument shows that each T2o-,^(xi, 0:2) decays in the same way, and in particular that it is 
bounded by C| (771772)'^ |- Using (3.17) we see that one can choose the constant C to grow 
with yU at most as some polynomial whose order depends on D through the power n in (3.17) 
required to show convergence of the mode sums. 

To summarize, the 0{g) correction to the time-ordered 2-point function may be written 

(T0^(xi)0„(x2))i^^ = -gfdkpG „{xi,X2) + T2^,<^(xi,X2). (3.20) 

The first term is precisely the 0{g) correction to the Hartle-Hawking state [3 i]. At sufficiently 
late times \rji^2\ 3> |?7lmax the second term decays like |?7i?72|'^; i.e., exponentially in the usual 
global time coordinate t. Furthermore, the coefficient of the decay term grows at large a no 
faster than a power law in a, and similarly for all T2a,^i{xiiX2) with Re < Re ex. 

3.1.2 The 0{g^) correction 

Let us now compute the 0{g'^) correction to the time-ordered 2-point function. Although 
it is somewhat tedious, the effort will be worthwhile as we will make use of this result in 
studying the 1-loop 0^ correction in the next section. The 0{g'^) correction is depicted in 
Fig. 3 (b) and is given by the expression 

{T(l)^^{Xi)(j)^^{x2))^^^ 



i / _9iy)9iy)\GaAy^xi)Ga2iy,X2)G^^{y,y) - Wa,iy,xi)G^^{y,X2)Wa3iy,y) 

+ W,, {y, xi)W^, {y, X2)Gl^ {y, y) - G,^ {y, xi)W^, {y, X2)W^, {y, y) 

(3.21) 
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This expression may be organized into two terms, each of which contains the 0{g) corrections 
to an appropriate 2-point function: 



(1) 



(3.22) 



The 0{g) correction to the time-ordered correlator (first hne of (3.22)) was studied in 
section 3.1.1 above, and the 0{g) correction to the Wightman correlator (second line of 
(3.22) can be analyzed similarly. In particular, using manipulations analogous to those that 
led to (3.10), one obtains 

(0<,,(a;i)0<,2(a;2))J^ =i [ g{y) \W^^{xi,y)G^.,{x2,y) - Gl^{y,xi)W^^{y,X2)\ 
1 



y 



M2 - M2 [f(^l)^'^2(^l'^2) - g{.X2)W„,{Xi,X2)] 



+ 



^'9{y)) K,(xi,y)V^a,(x2,i/) -G* (i/,xi)V^l^.,(t/,X2)l(3.23) 



After inserting (3.11) and (3.23) into (3.22) and rearranging terms one can again recognize 
the 0{g) corrections which we have already computed: 

(r0<^^(Xi)0<^2(^2))i^'' 

g{xi) {T(j) 

0-3 1^1 J 00-2 (2:2))*^ - (T0^,(Xi)0<^2(x2))i^ ^ 



+ 



(0^3 (2/)0a2 {.X2))^l^ ^^PVai (y, Xl] 



(3.24) 



Here (T0o-^(xi)0o-2(3;2))^ ^ denotes the same integral expression (3.4) as (T0o-i(a;i)0<T2(a;2))^^ 
but with a giy) replaced by g^{y)- Once again we may use (3.11) and (3.23) to simplify this 
expression. The result may be written 



(r0^,(xi)0^2(x2))J^ = gj 



Gg-l (Xl, X2) 



G^^i^i, X2) 



{Mf - M^){Mf - Ml) {Mi - Mf){Mi - Ml) 

G^^{Xi,X2) 



+ 



(Ml - Mf)iMl - Mi 



3(Tl(J2(T3 



(xi,X2), (3.25) 



with Tso-io-ao-al^i' ^2) the collection of integration terms 



T3aia2a3{xi,X2) = {gV g terms) + {VgVg terms). 



(3.26) 
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where 

{gVg terms) := 



1 



{V^g{y))g{x,) 



+ (VV(l/)) 
+ {^'9{y))9{y) 
+ {V^g{y))g{x2) 



(Mi - Ml) 

Ga,{y,xi)V'ij,Ga^{y,X2) - Wa,{y,xi)V' f,Wa^{y,X2) 
{Ml - Ml) 

Ga,{y,xi)V'f,G^2{y,X2) - Wa,{y,xi)V' f,W^^{y,X2) 
(Mi - Ml) 

Ga^iy,xi)V^ ^G^^{y,X2) - W^-^{y, xi)V' ^,W^^{y, X2) 



(Ml - Mi) 



(3.27) 



(VgVg terms) := 



{Mf-Mi){Ml-Mi)^^^^ 



W^<7i y) V^G<^3 (y, y)^vW„^ (x2, y) 

- W^^ixi^y)"^ ^W„^{y, y)'^^jW„^{y, X2) 

- W^^iy.xi)"^ ^W^^{y,y)'^,Ma2{x2,y) 
+ {y, xi) V^G* (y, y) V^l^„, (p, X2)\. 



(3.28) 



Let us simplify the lengthy expressions (3.27) and (3.28). The terms in (3.27) are of the 
same form as T2aia2{xiiX2) above; indeed, after a few simple manipulations we may write 
(3.27) as 



{gVg terms) 



{Ml - M|) ^'"^'"^ ^""^^ + (M| - M|) ^'"^'"^ '''^ 



(3.29) 



2(M| - Ml){Ml - Ml)^' "^"'^ 
Here 72/0-^0-2(2^1) 2:2) denotes the same integral expression (3.12) as ^2 0102(2^1 5 2:2) but with 
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g'^{y) in place of g{y)- The (VgVg terms) may be decomposed into the sum of the two terms 
Twa,a,,as{xi,X2) := - Ml) {Ml - Ml) 
Re yjj^V>'g{y)){V^g{y))W,,ix,,y)^,W^,{y,y)Vj,W^,{y,X2)^ , (3.30) 



(M2-M|)(M|-M|) 
Re ly y(V'^(7(y))(V'^^7(y))W^.,(^i,y)^MG.3(y,2/)^FVr..(^2,y)} . (3.31) 

Expanding the Green's functions in modes, one easily obtains 



''"''^bl/^" 1 + ^ - 2)x:,.3(L)x.3..(^)^..L(mXL(^2)Cf "'^/'(xi • a;2j 

.L=0 



(3.32) 



This leaves only Tgo-io-jo-s (2^1, 3^2)1 which one may treat similarly using G„{y,y) = 0{rj 
ri)Wa{y,y) + 9{f] — ri)W*{y,y) and the mode sum (2.11). One finds 

{Xl,X2) 



f2-D_ 

" 27r^/2(M|-M. 



L=0 



X / drjil + v')^''-'^/'g\m<Liv)d^K^MKL.Uv) + (^1 ^ ^2, 



(3.33) 



where 



Cl.,.:M = ^^2 _ ^2) j_ c^r^(l + r/^)(^-^)/V(r]) [<L(^)5n«.3L(r])J . (3.34) 

Note that the integral in (3.34) converges since the integrand is smooth and g'{ri) has compact 
support. Using (2.16), we see that the large L limit of (3.34) is a smooth function of f] which 
is in fact independent of L,Mi,M^. As a result, the f] integral in (3.33) once again gives a 
function of L which decreases faster than any power of L. In particular, as with Xo-io-2(-^)) 
for any p > it may be bounded by Cp{ai,a^)L~'P where Cp(cri, a^) is a polynomial in cxi, cxs 
whose order and coefficients are determined by p. 

It is clear that (3.32) and (3.33) are similar in form to (3.12). As a result, T^^ 0-1 0-2 0-3 (2^1, 2:2) 
and TGaicT2az (a^i, X2) can be shown to decay like rjl^rj"'^ via arguments analogous to those used 
for T2o-i(T2(a^i5 2:2). Furthermore, just as with T2aicr2{.xi,X2), we find that these functions are 
in fact bounded by C*| (771)^(771)2 1 for some polynomial function C(cri,(T2). 
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9? 



Figure 3: The 0{g'^{x)) correction to the time-ordered 2-point function (T0(xi)0(x2))^ in 
a theory with an g{x)(f)^{x) interaction. Once again we label each leg of the diagram by a 
distinct mass parameter aj. 

Collecting our results, we conclude that the 0{g'^) correction to the time-ordered 2-point 
function is 

G'ctj(Xi,X2) G'ct2(Xi,X2) Gfj^{Xi,X2) 

{Ml - Mi){Ml - M^) {M^ - Mf){Mi - Mi) {Mi - Ml){Mi - M^) 

^{Mf-Mif''"'''''^'''''''^ + {M^-Mif'''^'"-'^'''''''^ 
{Ml + Mi-2Mi) 

~^2(M2 - M'^){M'^ - M^) ^'"^"^ ^ +Twaia2,a3[Xi,X2) + Tq aia2,a3{Xi, X2) ■ 

(3.35) 

We remind the reader that T'2o-jo-2(xi, X2) is defined in (3.12), T2/ 0-^0-2 (xi,X2) is (3.12) with 
g'^{y) in place of g{y), T'vj/ 0-2,0-3 (a;i,X2) is defined in (3.30), and TGcTia2,cT3{xi,X2) is defined in 
(3.31). In (3.35), the term in square brackets is the 0{g'^) correction to the Hartle-Hawking 
state. The remaining terms are bounded for all xi, X2 G J^{TZ) /TZ and in particular are finite 
at coincidence {xi = X2). In addition, these terms all decay like rji^rj2^ when jr^i 2I > |?7|max- 
One may readily verify that the full expression (3.35) is regular in the limit of coincident 
masses. In this limit, and at late times, we obtain 

(T0o(xi)0.(x2))J^ = ^gjdl2G„{x^,X2) + O {{mmV) ■ (3.36) 

3.2 Example: (j)^{x) interaction 

In this section we consider a cubic self-interaction 

Ant[0] = -^0'(x), (3.37) 

and compute the 0{g'^) correction to the time-ordered 2-point function. The relevant Feyn- 
man diagram is shown in Fig. 3. In spacetime dimension D > A this correction contains 
ultraviolet divergences, so we will need to regulate our computation and renormalize the 
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theory. We restrict attention to -D < 6 for which 0^-theory is power-counting normahzable; 
for these dimensions, and to 0{g'^) the counterterms needed for renormahzation are 

= _(Mi^V,^(x)V>(x) - iM^l^Ml^^,), (3.38) 

with Z^{x) and Zm{x) given by Zi = 1+0(5'^). The renormahzation coefficients are position- 
dependant as a result of our position-dependant couphng g{x). However, no renormahzation 
is required for D = 2, 3. 

To regulate our computation of the diagram (Fig. 3) we replace the internal Green's 
functions with Pauli-Villars regulated Green's functions. Thus the full 0{g^) correction is 
given schematically by 

(T0,,(xi)0,,(x2))J^ = (diag) + (c.t.), (3.39) 
where (diag) is the regulated Feynman diagram given by the expression 

(diag) [ I ^7(yi)^?(y2)|G..(yi,a;i)G.,(l/2,X2)G^;3*^(yi,l/2)G^.7(2/i,y2) 

- W,, (yi, x{)G^,,{y2,X2)Wlf{y,, y2)Wlf{y,, y^) 

+ (2/1, x{)W„, (2/2, x2)GlY{yu y2)GlY{yi, 2/2) 

-a,(2/i,xi)l^.,(2/2,a;2)W^;:^(y2,l/i)W^^f(l/2,Z/i)},(3.40) 

and (c.t.) are the counterterms generated by (3.38): 

(c.t.) = -{Z^{X2) - l)G„^{Xi,X2) 

+tMl I {Z^{y) + ZM{y) - 2) [G,,(2/,xi)a,(2/,X2) - 1^^,(2/, xi)l^.,(2/, X2)] (3.41) 
•Jy 

At the end of our computation we will set the masses to be equal. 

A useful way to proceed is to make use of the linearization formulae for Bunch-Davies 
Green's functions: 

Hl\^{x,y)HlJ{x,y) = ! f{f,)H^{x,y). (3.42) 

Here Hcr{x,y) and Hl^^{x,y) may be taken to be any Bunch-Davies Green's function (e.g., 
time-ordered, Wightman, etc.). In addition to /x, the function /(/x) implicitly depends upon 
(Ji, (T2, and the spacetime dimension, as well as the collection of Pauli-Villars masses for 
the two Green's functions. To derive an expression for the function /(/i) it is sufficient 
to construct the linearization formula for the Euclidean Green's function Acr{x,y). Since 
the Bunch-Davies Green's functions are given by the Euclidean Green's function with an 
appropriately chosen prescription for avoiding the cut in the complex Z plane, the extension 
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of the linearization formula to these Green's functions follows immediately. We derive the 
linearization formulae in Appendix C following [39]; here we will simply state the results. 
For the time-ordered Green's function we have 

GlJ{x,y)GlJ{x,y) = e-"" [ (/i) + r-(/x)) y), P = 2,3) 

Gl\^{x,y)GlJ{x,y) = f-"" [ (/i) + G^(a;, y) - z£^-^Co5(a;, y), p = 4, 5) 

Gl\%x,y)GlJ{x,y) = f-"" [ (/i) + r-(/.)) y) 

-it-^Co5{x, y) - U^-^CiU5{x, y), {D = 6, 7). (3.43) 

where /o-icr2(/^) depends only on cri,cr2 and is independent of the regulator masses. In the 
complex n plane the function (/o-io-al/^) + Z™"^!/^)) decays exponentially away from the real 
axis and is analytic in the strip Re (di + (72) < Re /i. The contours of integration (3.43) and 
(3.44) lie within the strip Re (di + (T2) < Re /x < 0. 

The coefficients cq and ci in (3.43) are real functions of the Pauli-Villars masses but do 
not depend on x, y or the integration variable /z. These expressions have been organized to 
make it easy to take the limit of large Pauli-Villars regulator masses. As discussed in the 
Appendix, in this limit the function faia2{tA remains, /™(/i) vanishes, and Cq and Ci diverge. 
The explicit expressions for faicT2{f^)y /™°(/^)! co and ci can be found in the Appendix, but 
we will not need them here. For the Wightman Green's function we have 

W;i%x,y)W^,l^{x,y) = f^"^ [ {U.M + ^''{^^))W,{x,y), (all D). (3.44) 

This expression is finite when the Pauli-Villars masses are taken to infinity. This refiects that 
fact that the product of such Wightman functions (x, y) ■ ■ ■ (x, y) is positive /negative 
frequency with respect to x/y. 

Let us start with the simple D = 2,3 where there are no ultraviolet divergences at this 
order. In this case we may set Zm = = 1 + 0{g^) and immediately take the limit of 
large regulator masses^ (which sends /™'^(/^) to zero) in the linearization formulae (3.43) and 
(3.44). The full correction to the 2-point function (3.39) then becomes 

(T0.,(xi)</).,(X2))J^ = f-"" 

G„^{yi,xi)G„^{y2, X2)G'^(i/i, 1/2) - W^^ (yi, xi)G„^{y2, X2)Wf,{yi,y2) 



i / / 9iyi)9iy2 
yi Jy2 



+ W„^{yi,xi)W^^{y2,X2)Gl{yi,y2) - G^^{yi,xi)W^^{y2,X2)W^{y2,yi) 

(3.4 



^We could also have done the full computation without ever involving regulators. The unregulated Green's 
functions for Z? = 2, 3 satisfy the analogues of (3.43), (3.44) with = 0. 
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The astute reader will recognize the term in braces as the expression for the 0{g^) correction 
to the 2-point function in the 0^-theory discussed in §3.1. Inserting our final expression (3.35) 
for that correction we obtain 



(T0 (3^2))^^ 



02~D 



j /<T3^4(/^)|^/ 



Gfj^{Xi,X2) G(j^{Xi,X2) G^{Xi^X2) 



[{Ml - Ml){Ml - M^) {Mi - Ml){Mi - M^^) {M^^ - Mf){M^^ - M^ \ 
+ T3„^^2f,{xi,X2) > (3.46) 



— p'2~D [ f f \ \ 2 (^^(Xi, X2) ^2^(72(^1)^2) ^2ai^(a^l) ^2) 

J, r {Ml - Ml) {Ml - M|) + {Ml -MD^ (M| - M^) 

+ 7V<^,02M(a;i,a;2) + Tq a^a2Pi{.xii ^2)^ , (3.47) 

for x,y G J^{1Z)/71. To obtain the second line we have inserted the definition of 
(T2At(^i) ^2) and noted that terms in the integrand of (3.46) whose only dependence upon 
/i is a factor of l/{Ml — Ml) make no contribution to the integral over /x. This is because 
l/(Mf — Ml) contributes only poles to the left of the integration contour. For these terms 
the integration contour may be closed to the right. But there are no poles contained in the 
right half-plane, so these integrals vanish. 

In higher dimensions the only additional complication is that we must take care to cancel 
the ultraviolet divergences contained in (3.40) with our available counterterms. Our use of 
Pauli-Villars regularization as well as our linearization formulae make this procedure quite 
transparent. Consider first the case of D = 4, 5. After utilizing our linearization formulae 
(3.43) and (3.44) the divergent terms in {T(t)fj^{xi)(t)^^{x2))\, are 

+ cot-''i j g\y)[G„,{y,x^)G^,{y,X2)-W„,{y,x^)W„,{y,X2)]. (3.48) 

Comparing this expression to our counterterms (3.41) we see that these terms are cancelled 
by setting 



Cot'' 
Ml 



Zm{x) = 1 - ^^g\x) + 0{g^), Z^{x) = 1 + 0{g^), {D = 4, 5). (3.49) 



In D = 6 dimensions the divergent terms are 

- Cig)G„^{xi,X2) + {co£~^ + CiM^) i / g'^{y) [G^^{y,Xi)G^^{y,X2) - W^^{y,Xi)W^^{y,X2)] , 

Jy 

(3.50) 

which may be canceled by setting 
Zm(x) = 1- (^ + ci^/(x) + 0(/), Z^{x) = l-cig) + 0{g'), (D = 6). (3.51) 
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With all divergent terms cancelled by the counterterms we may take the Pauli-Villars masses 
to infinity, so that — )■ 0, and then proceed as for D = 2,3 above. For all 2 < D < 6 
we find 

{TMx^)M-2)h =i UM [9f _ + {Ml - M,^) + {Ml - M,^) 

(a;i,X2)|, (3.52) 

where we have now taken the limit in which all masses become equal. Note that the 
units are correct: mass^ has units g^{x^ has units and Gfj{x\,X'2), T2axa2{.^\-,^'i) 

Tw u^a2A^\^^2), TGaia2fi{xi,X2) each has units 

The first term in (3.52) gives precisely (T0o-(a;i)0o-(a;2))^^g^, the associated correction to 
the correlator in the Hartle-Hawking vacuum of the theory with g = gf = const, for all time. 
Each of the remaining terms inside the braces was analyzed in detail in section 3.1. Choosing 
the /i contour to satisfy Re/i < Re a, we see that each such term is bounded by C(/i)| 771772!°^. 
Recalling that C depends at most polynomially on fi while f^a^ decays exponentially at large 
imaginary /i we see that the integral of these terms over /i is bounded by C\r]i'r]2\^ for some 
constant C; i.e., 

2-D f f- f \ Gf^{xi, X2 



{Tct>^{x{)U^2))f = g)e-^ I /..(/i) (^"2 _'M2)2 + ^ ((^1^2)'^) , (3.53) 

where the first term is precisely the one-loop correction [■•'■] {T(j)^{xi)(j)(j{x2))'f^u to the 
Hartle-Hawking vacuum of the theory with constant g = gj. 

We emphasize that the O ((771772)'^) term is finite at coincidence {xi = X2). It follows 
immediately that the one point function of the composite operator 0^(xi) in our state can 
again be written as the sum of two terms, the first being its (finite) value in the Hartle- 
Hawking state for constant g and the other decaying like 'qf'^. In particular, defining (j)'^{xi) 
via any de Sitter-invariant regulator gives a result of the form 

(0'(^i)>l'^ = const. + 0((r^i)2'^). (3.54) 



4 Discussion 



We have shown that, in the time-dependent model of section 3, the two-point function at 
one-loop level approaches that of the constant coupling Hartle-Hawking state in the limit 
where its arguments are evaluated at late times 771,772- A key feature of this model is a 
time-dependent cubic coupling g{x)(j)^ which vanishes before some fixed initial (global) time 
which we may call t]q. We also chose g{x) to be a constant {gf) to the future some Cauchy 
surface. The state was taken to coincide with the (free) Bunch-Davies vacuum in the region 
where g{x) = 0, so that the state in the region with g = gf is determined by the particular 
way in which the coupling turns on. Since we required only that g{x) be smooth in this 
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transition region, this scenario describes a wide range of possible states. Our results were 
established by working entirely in Lorentz signature and in global coordinates. 

Our main result is that the difference between our two-point function and that of the 
g = gf Hartle-Hawking state is negligible for rii,ri2 ^ |?7o|- In other words, while the 
discrepancy may be significant (and perhaps even large!) for some period of time, it decreases 
rapidly once the universe enters its expanding phase and the size of the spheres becomes 
significantly larger than they were when the coupling turned on. This is precisely what one 
would expect based on the free theory in which perturbations rapidly disperse as, after this 
time, their effect on local correlators decays rapidly. This behavior was shown in [38, 39] to 
hold to all orders in perturbation theory for a dense set of states; our results here indicate 
that this dense domain allows for physically interesting initial conditions. 

The above model was recently considered by Krotov and Polyakov [!()]. While they 
characterized the model as "unstable," we remind the reader that their technical results 
are completely consistent with ours. As stated in their paper (below their equation (17)), 
their analysis applies in the regime rjo <^ — — 1^?2| (our notation); i.e., in precisely the 
complimentary regime to that studied here. As noted in the introduction, the divergence 
they find as r/o ^ ~oo is not only consistent with, but in fact is naturally expected from, 
the results found here and in [38, 39]. 

Despite various technical features of our analysis, we see that the approach to the Hartle- 
Hawking state at late times followed from a few simple ingredients. First, for quadratic 
perturbations g{x)(f)'^, one can write the full two-point function at each order as a sum of 
the Hartle-Hawking two-point function and a 'boundary term' associated with the transition 
region (see (3.9), (3.24)). These manipulations involve only integrations by parts and will 
clearly hold in a general spacetime which is asymptotically de Sitter to the future. The rapid 
expansion of the universe at late times then i) causes any given mode to decay as a power 
law in rj and ii) implies that the modes which have not yet decayed at some late time r) 
correspond to very high L. As a result, at the early time rjo when g{x) was time dependent, 
these high L modes were very high frequency. Since quadratic perturbations do not lead to 
loops, the Green's functions that appear in this boundary term are all positive frequency, at 
least at the key step (see the discussion of TGaia2,a3{xi,X2) surrounding (3.33)). As a result, 
in the boundary term these modes appear with coefficients involving what is effectively the 
Fourier transform of g{x) at large momentum, which vanishes rapidly since g{x) is smooth. 
Thus the effect of the boundary term decays with time, leaving only the Hartle-Hawking 
term in the two-point function. 

For the cubic perturbation g[x)(j)^, we used the linearization formulae (3.42) to make 
renormalization straightforward and to reduce the one-loop calculation to the quadratic- 
perturbation calculation described above. It is clear that analogous results follow immedi- 
ately in any context where similar linearization formulae hold for the associated free Green's 
functions. While such formulae are not obvious for general fields in general spacetimes, they 
must hold for conformally coupled free fields in general spherically symmetric spacetimes 
(which are necessarily conformal to dS), at least after inserting powers of the appropriate 
conformal factor Q{x). Note that this argument requires conformal invariance only for the 
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free theory about which we perturb'^ and that there is no restriction on the interaction. The 
only impact of these extra factors of fl{x) is to provide what is in effect an extra time- 
dependence in the couphng. Thus, to the extent that one can study a theory of general 
spacetime-dependent mass m^(a;) by perturbing the free conformally coupled theory by a 
quadratic perturbation g{x)(f){xy, our cubic results extend for all > to any spherically 
symmetric spacetime which is asymptotically de Sitter in the far future. 

It is useful to comment on the special case where the spacetime is taken to be the Einstein 
Static Universe (ESU) S*^ x M to the past of some S^. Since the ESU is static and spatially 
compact, we can take a limit where the coupling g{x) is turned on adiabatically slowly (after 
which it is = Of = const.) and in the distant past. We are then guaranteed that, in 
the ESU region, the state is given by the interacting ESU vacuum. As a result, subject to 
the same qualifiers as above we may consider the theory with g = Qf = const, for all time. 
Taking the state to be the (interacting) ESU vacuum at early times, we see that at late times 
correlators again approach those of the de Sitter Hartle-Hawking vacuum. This result can 
then be further generalized to either n— particle or thermal states in the ESU, all of which 
approach the same de Sitter Hartle-Hawking vacuum at late times. 

In summary, we have shown at the one-loop level that a wide class of physically-motivated 
initial conditions lead to two-point functions which approach that of the Hartle-Hawking 
state at late times. This suggests that states defined by general physical initial conditions 
lie in the dense set of states where the cosmic quantum no hair theorems of [38, 39] apply. 
We expect that this can be explicitly checked by extending the calculations reported here 
to all orders in perturbation theory. After all, the techniques used above were essentially 
Lorentz-signature versions of the Euclidean methods applied in [34, 39]. So by adapting 
further such techniques to Lorentz signature we expect to obtain all-orders results analogous 
to those found in [38, 39]. 

Some readers may feel a lingering uneasiness with these results due to the well-established 
infrared divergences (see e.g. [63, 26, 32, 30]) associated with in-out perturbation theory in 
global de Sitter. In particular, as noted in e.g. [ ;], at sufficient loop orders such divergences 
occur in the future expanding region even for very large masses. This certainly indicates 
that some quantity is becoming large in the infrared. However, the key point to realize is 
that the quantity need not be local. In particular, we suggest that it is merely the operator 
relating the (free) Bunch-Davies vacuum |0) to the interacting Hartle-Hawking state \HH) 
which becomes large at late times. This operator involves integrals over an entire at each 
time and can become large as the grows in size. In Minkowski space, a corresponding 
IR divergence is forbidden due to the exponential decay of massive propagators at large 
spacelike separations. But since the volume element also grows exponentially in dS, such 
IR divergences can occur. Indeed, the operator relating |0) and \HH) is closely related to 
the vacuum to n-particle amplitudes of in-out perturbation theory noted to diverge above. 
This stands in sharp contrast to the good IR behavior of (unintegrated) n-point functions 
as established here and in [34, 38, 39]. 



•^Recall that conformally coupled free fields correspond to cr = 



-{D- 2)/2 in de Sitter. 
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A Asymptotic expansion of Klein-Gordon modes 

One way to derive (2.14) is to use a standard Laurent expansion for the Gauss hypergeometric 
function. Recall that F„L{rj) is a Gauss hypergeometric function 



FaLiv) = 2F1 



D D ^ D 1-ir] 

a -\ ,1 — a \ L -\ ; 

2 2 2 2 



(Al) 



For 1(1 — iri)/2\ < 1 this hypergeometric function has the power series expansion: 



a + f,l-a-f J 



n=0 



o-+f + n, l-(T-f+n 
l + n,L + f+ n 
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while for |(1 — ir])/2\ > 1 it has the Laurent expansion: 
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(A.3) 



The asymptotic expansion of Fcriiv) fo^' large |?7| ^ 1 with a, D, and L fixed is given by 
(A.3); the leading terms are 
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(A.4) 



Sub-leading terms in the Laurant expansion (those with n > in (A.3)) are negligible when 
|?7| ^ 1 and \ri\ ^ (L — a). In this limit 



;i+^2)-(D-2)/4^^-(D-2)/2_^Q(^ 
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and so the full mode functions have the asymptotic expansion 



2cT+{D-2)/2 



L + f,2a + D- 1 
L + a + D-l,a+^ 



2 J 



exp 



1 + 



cr 



(a ^ -(ct + D- 1)), for > 1, jr^] > L - a. 



(A.6) 



B WKB approximation of Klein-Gordon modes 

In this appendix we derive the WKB approximations to the Klein-Gordon modes (2.16). 
The d'Alembertian in the coordinates (2.1) is 



(B.l) 



where V|_d_i is the scalar Laplacian on ^ with unit radius. Expanding the mode functions 
as (p^iix) = l^'^~^^/'^UaL{,'il)Yi{x) gives 



where primes denote derivatives with respect to rj and 



We now insert the ansatz 



u 



(B.2) 



(B.3) 



(B.4) 

and anticipate that A{ri) = Aq^t]) + Ai{ri) + . . . , with successive terms suppressed by a large 
parameter. Inserting (B.4) into (B.2) we obtain 

(1 + r/2) [A"ir]) + {A'iv))'] + Di^A{r^) + /(r/) = 0. (B.5) 

At lowest order in the WKB approximation we keep only the term (^0(77))^; this yields 



A'o(r^) = ±1 



;i + r^2 



(B.6) 



We see that the large parameter in this WKB expansion is \/ firj)^ and that we have just 
solved for the 0{f{ri)) part of (B.2). To obtain the sub-leading term Ai{r]) we solve the 
0(y7M) part of (B.2): 



1 + v') Kiv) + 2A',{v)A[{v)] + Dt^AM = 0, 



from which we obtain 



2 [1 + 7^2 A'^ir])] dr] 
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log [(l + r/)-^/^(A'o(r/))-V2] 
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(B.8) 
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We note that 



where^, 



^io(^) = (B.IO) 

is independent of L, M. This (and similar resuhs at higher order) shows that our WKB 
expansion is uniformly valid at large f{ri). 
Collecting our results we have 

u.L{r,) ^ ^(1 + ri'r^'"^/' [fm-'" e^^^(^), (B.ll) 

1 /2 

with T{r]) the anti-derivative of + rj^)] ■ The normalization is easily obtained by 

computing the Klein-Gordon norm (2.6) of (B.ll). The higher order terms give corrections 
to (B.ll) that are much smaller than (B.ll) in the limit of large /. 



C Linearization formulae 



In this appendix we construct linearization formulae for the Bunch-Davies 2-point functions; 
these formulae follow immediately from the linearization formula of the Euclidean Green's 
function. To construct the linearization formula for Euclidean Green's function we must 
recall some facts about the Euclidean Green's function ^„{x,y) = Ao-(Z), where Z := 
Z{x, y) is the SO{D + l)-invariant distance between x and y. This information is presented 
in more detail in [.39]. Recall that Afj{Z) may be written as a Mellin-Barnes integral 
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(4vr)"+i/2 

Here j^ - ■ ■ denotes a contour integral in the complex u plane traversed from —ioo to +ioo 
within the strip cr < Rez/ < 0. The measure dv/{2TTi) is assumed. Because the sphere radius 
C. enters only as a multiplicative constant we will set £ = 1 for now and restore it at the 
end of the appendix. It is convenient to use a := {D — l)/2 to keep track of the spacetime 
dimension. In Pauli-Villars regularization we subtract from the unregulated Green's function 
a linear combination of Green's functions organized such that all UV divergences cancel. We 
may write the regulated Green's function as 



(C.3) 



Here we used IF 



^1 

2/1 



< \F\ 



I /'I 

2L(L-D+2) 



(1 + 77^) for F{r]) 



Dri 



26 



with 



[D/2] 



(C.4) 



Here the are mass parameters corresponding to the PV masses M"^. The values of the 
coefficients Cj are not unique, though their values are constrained such that the desired 
cancellations occur, and they are bounded in the limit Mf — )■ oo [ ]. For finite PV masses 
the function iIj]^^{v) is analytic in z/ in the strip Re a < Re z/ < ^ in all dimensions. At large 
^ 1 the function ipai^^) has the asymptotic behavior 



M 



2a-l+2u 



(47r)"+i/2 
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a — u 
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(C.5) 



it follows that in the limit — )■ oo the regulated Green's function reduces to the unregulated 
Green's function at finite separations. Further details may be found in §3.1-2 of [39]. 
The product of two regulated Euclidean Green's functions is given by 



ji/i Ji/i 

One can easily invert (C.l) to find 



C?(^i)C^(^2)r[-z.i,-z.2] 



1 - z 



U1+U2 



(C.6) 



2(47r; 



0+1/2 / Y 



^ + a + u,fi — u 
1 + 2a + fi + u, — z/ 
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Combining (C.6) and (C.7) we immediately obtain the expression 



A-s(Z)A-/(Z) = / /(p)A,(Z) 



(C.8) 



with /(/i) given by 



/(/i) := 2(47r)"+i/2(/i+a) 



C?(^^i)C^(^2)r 



-Z/l, -Z/2, 2 + a + Z^l + Z^2, - '^l - ^^2 

-z/i - z/2, 1 + 2a + /i + z/i + z/2 

(C.9) 

In the complex plane the function /(/i) decays exponentially away from the real axis, and 
is analytic in the strip Re (cxi + (72) < Re/x. The contour of integration in (C.8) lies within 
the strip Re (ai + (T2) < Re /i < 0. 

Next we perform some simple manipulations, the utility of which will become clear later. 
Let us denote the Pauli-Villars parameters associated with ^^^HZ) by Ch, M^^, and p 



li 



respectively, and likewise those parameters associated with A^°^(Z) by 6*2^, M|j-, and p2j- 
Let us further define 



fp,pM := 2(47r)°+i/2(/x+a) 



^/3i(l^l)^/32(z^2)r 
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(C.IO)" 
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(this is just (C.9) with unregulated i'ai^) functions). By expanding the ip^^^iyi) as in (C.4) 
we may write 

[D/2\ [D/2] [D/2\ [D/2\ 

i=l j=l i=l j=l 

Consider the terms fp^^p2j (/^) iii the last sum of this equation. In these terms we would like 
to move the contours of integration into the region Re z/i < | — a, Re 1^2 < | ~ For D > A 
we can only do so at the cost of encountering poles. Performing this manipulation yields 
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2(/i + a). (C.12) 



Here hp^^p^.i^^) is just /piiP2j(/^) with the contours satisfying Rei^i < ^ — a, Re 1^2 < ^ — a as 
desired, and the sum is due to the residues 
are coefficients that do not depend upon /i: 



desired, and the sum is due to the residues of the poles encountered for D > A. The c" 
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Defining yet another quantity. 



[D/2] [D/2] [D/2] [D/2] 

r-(/i) := J2 Cufp^.M + E c2,u.pM + E E Ci.c^2,v.p.,(/i), (C.14) 

i=l j=l i=l j=l 

we may write /(/i) as 
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None of these manipulations have altered any of the salient features of /(/i). 

The purpose of these manipulations has been to isolate the part of /(/i) which diverges 
in the limit where the Pauli-Villars regulator masses M?j — )■ 00. In this limit: 

1. the function /o-i(T2(/^)) which is independent of the regulator masses, survives unaltered, 

2. every term in decays like a negative power of at least one regulator mass, and 
so it vanishes in the limit, 

3. the terms in the sum on the last line of (C.15) diverge like positive power of a regulator 
mass. 
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In particular, for the spacetime dimensions of interest for our computation we may write 
these divergences exphcitly: 



/(/^) = /....(/i) + /n/^), 

= U^M + /™"(/^) + Co2(/i + a), 
/(/w) = fa^M + + co2(/i + a) + ci2(/i + a)[-fi{fi + 2a)], 



P = 2,3) 
(D = 4,5) 
(D = 6). (C.16) 



Here cq and ci are coefficients that diverge with a Pauh-Villars mass hke cq ~ logM in 
D = 4, Co ~ M in D = 5, and cq ~ logM and ci ~ logM in D = 6. We have not yet 
taken the M^j — )■ oo hmit, but we have organized /(//) so that this hmit will be quite easy. 
We are almost done. Noting the relations 
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and restoring the radius i, we may finally record the linearization formulas for the dimensions 
of interest: 
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The linearization formulae for Lorentzian Green's functions may be found by the usual 
analytic continuation. For the time-ordered Green's functions we analytically continue 
Z ^ Z + ie. The resulting linearization formulae are given in (3.43). To obtain the lin- 
earization formulae for the Wightman Green's function we analytically continue with the cut 
prescription Z ^ Z + s{xi,X2)ie, where s{xi,X2) = +(— ) if xi is in the future (past) of X2. 
This cut prescription removes contact terms, so the result is the simple expression (3.44). 
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